Nonlinear dielectric susceptibilities in supercooled liquids : a toy model. 
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The dielectric response of supercooled liquids is phenomenologically modeled by a set of Asymme- 

tric Double Wells (ADW), where each ADW contains a dynamical heterogeneity of Ncorr molecules. 

We find that the linear macroscopic susceptibility xi does not depend on Ncorr contrary to all higher 

order susceptibilities X2fc+i- We show that X2fc+i is proportional to the fc*'^ moment of Ncorr, which 

could pave the way for new experiments on glass transition. In particular, as predicted by Bouchaud 

■ ^ ' and Biroli on general grounds [Phys. Rev. B, 72, 064204 (2005)], we find that xs is proportional to 

^ , the average value of Ncorr- We fuUy calculate X3 and, with plausible values of few parameters our 

^— ^ . model accounts for the salient features of the experimental behavior of X3 of supercooled glycerol. 

(n: 

H \ Upon fast enough cooling, most liquids do not cristallize but enter into a supercooled liquid state [Íl-[al, where the 

I I I viscosity 77 dramatically increases with lowering the temperature T. Below the glass transition temperature Tg, r¡ ia 
^. ' so high that the system is in practice a sohd -the glass-, yet, no structural difference between the glass and the liquid 
state has ever been detected [l|. Over the past fifteen years, a major breakthrough was the discovery of Dynamical 
Heterogeneities (D.H.) in supercooled hquids |6l-[ll| ; i.e., relaxation happens through coUective events gathering Ncorr 
I— 1¡ molecules, and some groups are relaxing much faster than others. As it is expected that an increase of Ncorr when 
^^ ' lowering T could increase dramatically 77, an significant efFort was made for measuring the T-dependence of Ncorr 

P: BülQ. 

It has been argued that the most direct way to draw accurately the T-dependence of N^orr from experimental data 

is based on the a.c. nonhnear susceptibility xa ¡14| - [l6| . where X3 is the third order response of the fluid to a field 

with an angular frequency lü. This field can be of any nature, e.g. electric as in 15, 16]. More precisely two nonhnear 

I \ susceptibilities are related to X3 ■ X3 ^^^ X3 j which correspond to the third order nonhnear response at the third 

' ^ ■ harmonics (i.e., at "3w") and at the first harmonics (i.e., at "Iw") respectively. Bouchaud and Biroh (BB) have shown 

Q . [3 [3 that X3 and X3 should be related to the average value of N^orr over the various D.H.'s existing at a given 

O ■ T -noted [N corr{T)\av- by : 
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fT^ \ Here k^ is the Bohznian constant, and Ta{T) is the typical relaxation time at temperature T corresponding to the 
^^ ■ relaxation frequency fa — 1/{2t:To,) where the imaginary part of the linear response is maximum. Axi = Xi{^ = 
^^ \ 0) — xi(i^ — >■ oo) is the part of the static linear susceptibihty corresponding to the slow relaxation process of interest, 
'^ . a^ is the volume occupied by one molecule, and "H and /C are two complex scahng functions that approach zero for 
¡> both small and large wTq. Note that the humped shapes of [H {uJTa) [ and [JC {u!Ta) [ are distinctive features of the 
glassy correlations. 

BB's prediction relies on very general grounds, such as a generahsed fluctuation dissipation relation, and was 
j3 \ inspired by spin glass physics [l7!|, where a true second order phase transition happens at Tc, accompanied by a 
critical divergence of X3 (while the linear suceptibihty xi does not divergc). A consequence of this generahty is that 
the detailed expressions of the scaling functions ?^ and K. remain unkown. Here we present a phenomenological "toy" 
model where, for the first time, BB's predictions are recovered with an explicit expression for the functions H and 
/C. By using plausible values of free parameters, the most salient experimental features of [l^, [l^ can be accounted 
for. Moreover we obtain new predictions on higher order nonlinear susceptibilities X2fe+i>5- This could motivate new 
experiments deepening our understanding of the glass transition. 

Model : We assume that aU D.H.'s are independent from each other and that a given D.H. is a group of N^orr 
molecules evolving in an Asymmetric Double Wefl potential (ADW) , depicted in Fig. [1] Each ADW is characterised 
by the height of its barrier V and by an asymmetry energy A. We neglect internal field effects. On Fig. [TJ z represents 
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FlGURE 1: Inset : ADW model where each D.H. 
with respect to the apphed ñeld E. Mam graph 



of Ncorr molecules evolves in an asymmetric double well making an angle Oi 
distribution of relaxation times for Glycerol at T = 204. 7K, lla. 1231]. 



the axis of the external electric field E{t) = E cos(a;í), and 6i is the angle between the field and the well which has the 
deepest energy at i? = 0. For simphcity we assume that O^ = 9i + tt. With respect to earher versions [l^ [20], a key 
refinement is the assumption that the magnitude of the net dipolar moment fi, in either of the two wells, is given by 
A* = ¡J-moiecy/ Ncorr whcrc firnoiec is thc molccular moment : This estimator of /i is assumed here because there should not 
exist any geometrical ordering among the molecules contributing to a given D.H. [6]. With vdh — NcorrO'^ the volume 
of a D.H., the simplest approach, for é'i = and A = 0, yields a static polarisation given by (/í/vd//) tanh{^E /kBT). 
Expanding in E gives Axi oc fi^ /vdh, which is independent of Ncorr since the Ncorr dependence of /j,^ cancels that of 
vdh- For ah higher orders such a cancenation does not happen, e.g. X3 °^ tJ-'^ /vdh oc Ncorr- This is the main reason 
why we find below that xi is bhnd to N^orr contrarily to all higher order susceptibihties. 

Let us now consider a set of N identical ADW's. With Hj^fc the transition rate from the weh k to the well j, 
we obtain the number m -resp. n^- of ADW's in state 1 -resp. state 2- by solving the two master equations : 
dni/dt = — H^irii -|- Hi^n^ and dn^/dt = —Hi^n^ + H^irii. Assuming thermahy activated barrier hoppings, one gets 
[¡113 : ni2 = VKexp[(A/2 + /íí;cos6'i)/fcBT], H21 = W^exp[-(A/2 + /i£;cos6li)/fcsT] where W = iyoexp[-V/kBT]. 
Here vq = 1/Tm where Tm is the microscopic characteristic time of the thermal fiuctuations within each wefi. The 
polarisation P of the set of N identical ADW's is given by P = /í cos(6'i)(ni — n^^/^NvDn)- The two niaster equations 
yield the dynamical equation for P, which involves the relaxation time r = 2W cosh{A/2kBT) of the identical ADW's : 
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P ((5 sinh e + cosh e) = A^ (ácoshe + sinhe) 



where e{t) = F cos{ujt) , F 
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Setting £; = in Eq. [5]yields P = P^ = AA6. As A^ ^ cos^i, we obtain < Pq >= where the brackets denote the 
average over the isotropicaUy distributed vahies of 6*1. In the hmit of smah fields (i.e., e — >■ 0), expanding Eq. [2]to the 
first order in e yields : 



„,, <MF>{l-5^) 
< Pi{t) >= — cos(a;í — arctanwr), 

^/1 + (wt)2 



(3) 



i.e., a Debye response, as expected in any double well model [21|. As the hnear dielectric spectra of supercooled 
hquids are asymmetric in frequency, we assume, as in other phenomenological models [llj, |22|, that the values of r are 
distributed according to G{t), [23]. Q{t), given in [l8[ and in Fig.[Tl is chosen to recover accurately the experimentaUy 
well known hnear susceptibihty xi (w, T) by weighting Eq.[3]with G{t) and summing over ah values of r. More precisely, 
G{t) determines the shape of xi(a;,T), but not its overall magnitude A^i. We use the experimentally weU known 
value of Axi as an additional constraint in our model : from Eq. [3l we obtain fjí^c.iec ~ SfcsTeoo'^Axi/^l — 5"^) ; i.e.. 
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FlGURE 2: (Color Online) Effect of the dimensionless asymmetry energy 5 on \X¡^ „J : the spectra have a low pass character 
excepted close to 5* — 1/y/Z. Here [Naorr\av = 5. 



^J'moiec is no longer a free parameter. 

Computing Xs ond x^ '■ For a given value of Ncorr and of S, we consider all ADW's having the same r and 6i. 
By using Eq. [51 we compute the polarisation up to the third order in field. We first average the result over ^i, then 
sum over r with weight G{t), and finally average over the values of Ncorr existing among various D.H.'s. The latter 
average is denoted by []av This yields quantities -labelled below by an index "m" standing for "model"- which are 
comparable to experiments [Ía lla]- Note that this method, by using the values of a/', Axi, Qij) drawn from standard 

(3) 

experiments, ehminates 9i,vo, V and fimoiec. Thus, when comparing our model to the experimental values of x\ £^nd 
Xg at a given T , the two remaining free parameters are \Ncorr\av and i5. For simpUcity we take a single value for 6, 
and postpone the possible averaging over 5 to Ref . [l^l • 

In practice, we solve Eq. [2] by assuming e <C 1, and develop Eq. [21 as well as P, in series of e, up to the third 
order. As the polarisation of a given set of ADW's sharing the same 9i,t, N^orr is not symmetric with respect to field 
reversal E — ^ —E, we set P{t) = X]o=o ^qi^) where Pq oc E"^. Since e ^ 1, one has \Pq\ » |íg'><j| ; i-e., all Pq'>q can 
be neglected when looking for Pq. Thus, Pq is obtained by keeping only the terms oc e' in Eq. [2l This was iUustrated 
above to get first Pq and then Pi{t) -see Eq. [31-. Repeating the procedure to the order e^ yields P2, [Í3|. Finally going 
to the order e^ gives : 



.d{P3l 
dt 



P^ 



X(l-á2) 

6 ' 



/■' 



SP^e. 



(4) 



As Pi and P2 are known, the analytical expression of Ps^t) is readily obtained from Eq. [H After averaging over 
9i, T, Ncorr-, onc obtaius Ps^m that must be identified with the third order term ^3 of the experimental polarisation. 
As e^ (X E^ {3 / 4: cos{ujt) + l/4cos(3ají)), we recall that ^3 naturally defines the first and third harmonics cubic 

susceptibiUties (with phases —¿3 , fc = 1,3) as [2J| : 
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E^ 



I cos(c.í - ¿W) -I- —\xf\cos{3u:t - S^/). 



(5) 



Results of the ADW model : A dimensional analysis shows that our model yields Pq^m oc [< M.F'^ >]av, which 
has two important consequences. First it yields Pq^m oc< (cos^i)*"'"^ >. This imphes that the even terms P^k^m = 0, 
which ensures that the macroscopic polarisation reverses exactly upon the E{t) -^ —E{t) reversal, as required by 
macroscopic symmetry considerations. Second, all odd terms P^fc+i.m are non zero, yielding for the susceptibihties : 
X2k+i,m oí [A"¿,^^] . This shows that the linear susceptibihty Xi.m is bhnd to the value of Ncom contrary to higher 
order susceptibihties which are directly proportionnal to the fc*'' moment of Ncorr- This first important resuh is 
reminiscent of the spin-glass transition [l7[ which has inspired BB's prediction. 



The above mentionned analysis yields x\ m 

3,miexp[ 



xflkBT/\e^a^{Axi)\ Writing X^^^ = 1^'^'^ 



that we convert into its dimensionless form X3 ^ = 
-iS\, ll where «^ 



-1, we get finally [í^l, with x 




FlGURE 3: (Color Online) For [Ncorr]av = 5, 5 = 0.60 comparison of the ADW model with the experiments of Ref. [I^ at 



T = 204. 7K. Xg ¿j is the weighted sum (see text) of Xg ^ (see Eq. |6]) and of Xg ¿¿„¿j,¡ corresponding to the cubic susceptibihty 

X3%t and the 



For f/fc > 1 one has Xf^ 



of independant molecules undergoing rotational brownian motion [161 . [2 

experiments are very well accounted for by the model. For f/fa < 1, only the global trends of the data are restored by the 

model with fex/fa = 0.14. Inset : Phases corresponding to the main graph, same symbols 
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FlGURE 4: (Color Onhne) Same symbols as in Fig. [3l excepted that Xg is displayed here, and that [Ncorr]av = 15. <5 and fs, 
have the same values as in Fig. |3l 
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Note that G nearly obeys Time-Temperature Superposition (TTS) ; i.e., it is nearly ¡22¡ independent on T when 

plotted as a function of A = r/Ta. As wr = XujTa, Eq. |ü] shows that X^ ^ equals [Ncorr]av times a function which 
does not depend on T -we take 5 as a constant in T-, when plotted as a function of wTq,. Thus Eq. |ü]gives the first 
phenomenological expression of the function T-LiuJTa) of Eq. [Tj-we recall that according to BB's prediction X)^ ' is 
[N corr]av'H- ■ Eq.[B]thus shows exphcitly that the T dependence of X^ is directly that of [Ncorr]avi up to smah effects 



coming from smah violations of TTS in Q{t). This is the second important result of our modeL 

Fig.[2]shows the frequency behavior of 1X3 ,„|. For most values of 5, the spectrum has a low pass character. In the 



vicinity of 5* = l/-\/3 the spectrum has a humped shape. To understand this, let us note p**"* the solution of Eq. [5] at 
w = 0. One gets P"*"* = M tanh[e + A/(2fcsT)]. Expansion to order e^ yields X|*^*(<5*) = 0. Around S*, X|*^* moves 
from a negative "Ising-hke" value (low á's), to a positive value for very asymmetric ADW's (high (5's). When w ^ 0, 
the effective relaxation time T/(coshe + ¿sinhe) comes into play, which contributes also to Ag ^. This is why close 

to S* , \X^^\ has a humped shape in frequency. A deeper, i.e., much less model dependent, reason for this humped 
shape is given below. 

To compare our model to the nonhnear susceptibihties of glycerol reported in [15|, [ig] , we first focus on the case 
/ > fa- Fig. [3]shows that choosing [Ncorr]av = 5 and 5 = 0.60 yields a very good agreement between our model 
and the vahies of X^ '{f > /„) measured at T = 204. 7K ~ Tg + 16K. We emphasize that the agreement is good for 

both the modulus and the phase of X^ ^. Fig. Ujshows the same kind of comparison for X^ , for which an expression 
similar to Eq. |6]is given in [l^. On Fig. |4]the best agreement between our model and the data reported in [lg| is 
obtained with [Ncorr]av = 15 and S = 0.60 : With respect to the data, our model underestimates the phase by ~ 20° 
and yields a maximum for the modulus at /* ~ 1.6/ti not far from the experimental value of 2.5/q,. The fact that the 
optimal [Ncorr]av is uot the same in Fig.Oand in Fig.|3]may come from interferences between the nonhnear responses 

of the D.H.'s with different r, see Eq. [6]and Ref. [l^. These interferences have different effects on X^ ^ and on Xg ^, 
see [l^, and this is not fully captured by our toy model, due to its simphcity. We emphasize, on the other hand, that 
[Ncorr]av = 5 — 15 is thc right order of magnitude when comparing to the values given by 4D-NMR experiments [^ 
or by Ref. [ÍJ|. Moreover our model accounts for the fact that \X^ \ is peaked at a frequency ten times smaUer than 
jXg I in glyceroL Finahy, S = 0.60 amounts to A ~ 1.4A;bT « lAkBTg, i.e. it does not introduce a new energy scale. 
Now, let us move to the case / < /„. Here we must take into account the ñnite hfetime t^x of D.H.'s; i.e., the 
fact that the hquid ñows at large times [T^]. The effective value of [Ncorr]av decreases with frequency when fxcx < 1, 
since a given molecule is involved in various DH's at large times ; i.e., it becomes independent of other molecules 
in the long run : as X^ ~ ' oc [Ncorr]av, this will give a humped shape to the nonlinear susceptibility even for the 
values of S where jA'g,^] has a low pass character. To take this idea into account, we simply use the -weh known- 
nonhnear response X¡n'^^^^^i of independent molecules (see [ÍalÍfllS^) and assume that it dominates the measured 
X^ when fTcx ^ 1. In practice, we write heuristicahy the total cubic susceptibihty X^ /^j (with once again A: = 1, 3) 



as 



^3%t = PX3% + (1 - P)^¡%tv^ai with p = exp i-fcx/f), see fl§¡. For / > /„, X¡%f is of course very close to 

3 „, since p ~ 1. For / < /„ Figs. Hashow that, with fcx = 0.14/„, X¡% 

(k 

3 



ATg l^, since p ~ 1. For / < /„, Figs. |5]13]show that, with fcx = 0.14/^, Xg ¿¿ has the same global quahtative trends 
as the measured X^ in glyceroL We note that fcx/fa = 0.14 amounts to Q = Tcx/tqi — 7, which is compatible with 
the values Q ~ 3 — 10 reported before [Ifl] albeit stiU debated [T^. We think that the osciUation of \X\ ¿J around 
0.1 /ct is unphysical and comes from the very naive way of inchiding t^x in our analysis. 

To conclude, we have developped a very simple toy model for the nonhnear susceptibihties in supercooled hquids. 
We find that X2fc+i,rn c>c [N^c.rr\a-v ; i-6-, that Xi,m is bhnd to the value of N^orr contrary to aU higher order suscepti- 
bihties. This yields the first phenomenological expression of the scahng functions involved in BB's predictions. With 
reasonnable values of parameters, the main trends of nonhnear experimental data are recovered. Our model explains 
very simply why the nonhnear responses yicld brand new information on the glassy dynamics. This simphcity may 
trigger more experiments deepening our understanding of the glass transition. 

We thank R. Tourbot for his outstanding technical help, S. Nakamae for carefuUy reading the paper. We thank G. 
Diezemann for encouraging discussions in dec. 2011, and C.Alba-Simionesco, G. Biroh, J.-P. Bouchaud, J.-P. Carton, 
P.M. Déjardin for long lasting help. 

Note added : See also on the same subject the paper of Gregor Diezemann to appear on Condmatt 
TODAY plus or minus a few days !... 
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We detail hereafter the calculations for Xg ¿^ and Xg ,„ summarized in our main article. We thcn give a short 
justification of our assumptionp = exp[— /gj;//] made in the end of the main article. Finally, we give more informations 
about what happens when averaging over the dimensionless asymmetry parameter 5. 

I. THE NONLINEAR SUSCEPTIBILITIES IN THE ASYMMETRIC DOUBLE WELL POTENTIAL 

MODEL. 

A. Calculations for one set of identical Asymmetric Double Wells. 

In this section we consider a set of N identical Asymmetric Double Wehs (ADW) ; i.e., a set of ADW's sharing the 
same valucs for aU microscopic parameters of the model. Denoting ni (respectively n^) the number of ADW's in state 
9i (respectively 62 = 01+71), the polarisation P of the considered set of ADW's is given by : 

(ni - n2)^imolecy'NcorrCOS01 n fXmolec COS 6i 

P == ^ ^TTTT 5 = M.— where M = — , (7) 

NNoorra^ N VK^a^ 

where it was assumcd that the net dipolar moment in either of the two states of a given ADW is given by /lí = 
^J'moiecVNcorr, with Hmoiec thc molecular dipole moment (see the main article). Combining the two master equations 
for ni(t) and n^^t), with n^it) = N ^ ni{t) , one gets the equation for the dynamics of P : 



dP 

T h P ((5 sinli e + cosh e) — A^ (ácoshe + sinhe) 

dt 



T. í,\ T^ r ,\ -n f^molecV ^corr COs(é'i) 

where e(í) = Fcos(wí) , F = ^ — -E 

kbT 

, .. t^molccCOsOi A 

ancl M = — ^ , ¿ = tanh(— — ). (8) 

As explained in the article, the two sources of nonhnearity in Eq.[5]are : (i) the nonhnear character of the equihbrium 
value P^*"* = A^(ácoshe + sinhe)/(ásinhe + coshe) = A^tanh[e + A/{2kBT)] ; and (ii) the nonhnear character of 
the instantaneous relaxation time Teff = t/{S sinhe + coshe). 

We expand P{t) in series of powers of the field E up to third order P{t) = Pq + Pi{t) + P^^t) + P^^t) where Pq oc E'¡. 
As E{t)^ = £'2(1 + cos(wí))/2 and E{t)^ = E^{3cos{ojt) + cos(3wí))/4, P^^t) and P3 are the sum of two terms : 

P2(í)=pf+pf(í) 
P3{t)=P¡'\t)+P¡'\t), (9) 

where the superscript in parentheses indicates the index of the relevant harmonics. For example, Ps^t) is given by a 
term osciUating at the fundamental frequency, and by a term osciUating at three times the fundamental frequency. 

As the condition e ^ 1 is weh obeyed experimentahy, one can neglect aU Pq'>q terms when computing Pq. Therefore 
Pq is obtained by keeping only the terms oc e^ in Eq. [5]above. 

To the order e° it is found that : 

Po = MS. (10) 

Now, going to the order e^, one has (by using the resuh for Pq in Eg. fTÜ)) : 

dP 

T — - + Pi=M\l-S^]xFcos{ujt), (11) 

dt 



which yields 



Md-S"^) 
Pi (í) = ^ ^ Frns {ut - arctan(a;r)) . (12) 

yi + (ZJ7)2 

We now go to the order e^ and get : 

T^^ + P2 = -SFPi{t) cos{Lüt). (13) 

As Pi{t) oscihates at frequency w, the right hand side of Eq. 1131 contains one constant term and anothcr tcrm 
osciUating at 2uj. Therefore, one finds : 

Pf = :^l^Í=£l£:cos[. + arctan(..)] 

P2^\t) = — — cos [2ijjt + TT - arctan(wT) - arctan(2wT)] . (14) 

2^1 + (cjt)2^1 + (2a;T)2 

FinaUy, we reach the order e^ and get : 



r^ + P3^ (1/6)M(1 - S')[e{t)f P2{t)5e{t) P,{t)[e{t)f /2. 



(15) 



(2)| 



We separate the terms oscillating at uj from those osciUating at Sw. Denoting by |Pi| -respectively [P^ \- the 
amphtude of Pi(í) -respectively P^ {t)-, one obtains : 



.d{pih 

dt 



Pl 



(1) 



(1/8)7V/ÍF^(1 - 5^) cos(wí) - (1/4)F2|Pi| cos[wí - arctan(a;r)] 
+ (1/2) cos[a;í + arctan(wT)] - PÍfhP cos(a;í) - 



>(2) 



(1/2)1^2 \5F cos[üjt + TT — arctan(a;r) — arctan(2a;r) 



(16) 



as well as 



AP¡'^) , .(3) 



dt 



+ PÍ 



{Í/24:)MF^{1 - 6^) cos{3ujt) - (l/8)p2|Pi| cos[3a;í - arctan(a;r)] 
-(l/2)|P2'^VPcos[3a;í + tt - arctan(a;r) - arctan(2a;r)]. 



By using Eqs. ll2lfT4l and the two previous equations, one finds 



(17) 



AP¡'') 

dt 



Pi 



(1) 



{IIA)MF^{1 - 5^)Df\ujT) cos Lí + «-^^^(a^r) 



where i:'^^(a;r)cos(a;í + í'_^^(a;r)) = - 



2^2 



1 

2^1 + (a;r)2 

cos [ujt + arctan(a;r)] 



cos(a;í) — 



cos [ut — arctan(a;r)] 



2yrTW 

5^ cos [ujt + TT — arctan(a;r) — arctan(2a;r)] 

^1 + {UJT)^^T+J2^ ' 



(18) 



as weh as 



APr) 



(3)^ 



dí 



^^+pf = (l/4)XP='(l-¿2)7^f (..r), 



l(3) 



^a^í + í-f^^a^r) 



with £)^''^(a;r)cos(3a;í + í'^''^(a;r)) = (1/6) cos(3a;í) 



(3)/ 



cos [3a;í — arctan(a;r)] 
2^1 + (^t)2 



(5^ cos [3a;í + tt — arctan(a;r) — arctan(2a;r)] 
^1 + {ujt)^^+{2ujt)^ ' 

(3) (3) 

Note that the above definitions of Dg and í'g arc consistent with thosc in the main article. 
The solution of Eq. [18] is given by 



(19) 



P^/\t) = "^^, ! ^'^^^^I d^^ujt) cos Lí + «-^'^(ojr) - arctan(a;r) 



4^1 + (^t)2 
The solution of Eq. [19] is given by 



Pf V) = — /^ '^ P>f\oJT) cos [3a;¿ + í-f ^(a;r) - arctan(3a;r) 

4^/1 + (3a;r)2 L 



(20) 



(21) 



B. Averaging over 6i,t and Ncorr- 

As explained in the inain article, three kinds of averages must be done in our ADWP model : 
(i) First, we have to average over the angle ^i, the values of which are assumed to be isotropicahy distributed. 
Denoting this average by < >, one finds < [cosOiY^^'^ >= and < [cosOiY^ >= l/(2fc + 1), for any integer k. As we 
have found above that Pq cx MF'^, one obtains < Pq >oc< A^F^ >cx< (cos^i)'^+^ >. Therefore, ah the even integer 
harmonics vanish, contrarily to ah odd harmonics which are found to be 



< P2fc+i >cx N¡:^„. 



(22) 



(ii) Second we have to average over various relaxation times r, with weight QÍT^dr. The distribution function Q is 
chosen so as to recover accurately the experimental hnear response, Xi- Therefore Q must simuhaneously solve the 
two following equations for the real part, xi, and the imaginary part, x'í ■ 



XÍ(w)-x'i(oo) 



Axi 



a(lnr) X 



1 + (wr)2 



Xi'(^) 



a(hir) X 



X dlTi: 



X dlnr. 



Axi y_^"^"'-^-- l + (cür)2 

where we have used the fact that Q{T)dT = ¿J(lnr)(ilnr. In practice one uses [22,123] : 



(23) 



Q{Ht)) = NcGEe 



with Ngge = 






TO 



l-P' 



(7-/3)/a_ 



(f)+-'-"(f)" " r(J) 



(24) 



Here T{x) is the Euler gamma function and a, /3,a,j,To are T dependent parameters. For glycerol, a good set of 
parameters is given by : 



a = 10 

f3 = -5.5996 X 10-1 +4.0900 X 10--''r+ 1.50795 X IQ-^T^ 



a = 1.57 X 10^1 exp 



407.525 



T- 141 



7 = -7.826920 + 1.015 x lO^^T - 4.32345 x 10"'' x T^ + 6.34415 x 10"^ x T^ 



ro = 1.1511 X 10""'" X exp 



19.08905 X 127.38588 



T- 127.38588 



(25) 



with T expressed in Kelvins. 

Note that ro is nearly proportionnal to the typical relaxation time t^ defined by t^ = l/{2T:fa) where fa is the 
frequency of the peak of x'i- Additionally one finds from Eq. [23] : 



Axi = 



(l-¿')(A^^oiec)'<[cOs(^l)]2> 

ksTeoa^ 



(26) 



and with < (cos^i)^ >= 1/3, we obtain 



yl^molec ) 



2 _ ^fcsTeoa^Axi 



(1-52) 



(27) 
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At this point the two free parameters of our ADWP model are Ncorr and 5. In particular Eq. [57] sets the value of 
^moiecj since Axi and a^ are experimentally well known. 

(ii) The third and the last average to be taken is over the vahies of N^orr- Indeed the proportionality expressed by 
Eq. [22]remains true when averaging over the r's. Therefore, it is very easy in our model to take into account the fact 
that there exists a distribution of the values of N^orr among various dynamical heterogeneities of a real supercooled 
hquid. As all above equations have been derived for given free parameters Ncorr and S, we superpose the ensemble of 
models with the same S but with different vahies of N^orr- Denoting the average over N^orr by [ ]at,, we obtain from 
Eq.HH: 



[< Í2fc+1 >]av OC X2k+l.,,n OC [iV* 



(28) 



where X2fc+i,m is the macroscopic nonhnear susceptibility of the order 2fc+l, and where the index m stands for "model", 
so as to avoid any confusion between the nonhnear susceptibihties produced by the model and those corresponding 
to what is experimentaUy measured (denoted x^fc+i)- Note that X2fc+i,m represents genericaüy the set of components 
of the macroscopic polarisation which is proportionnal to E'^^'^^ and oscihates at one of the odd harmonics between 
Iw and (2fc + l)aj. For example, X3,m corresponds to two terms : one is proportionnal to X3 m -note the presence of 
exponent (1)- and osciUatcs at Iw, and the other one is proportionnal to X3 m and osciUates at 3cli. 



C. Explicit expressions for the cubic susceptibilities Xs 1X3 • 



The macroscopic polarisation 3^ is given by [2J| : 



eo 



Xi{t - t')E{t')dt' 



X3(í - í'i,í - íí,,í - t'^)E{t'^)E{t'^)E{t'^)dt'^dt'^dt'^ + 



(29) 



where the function Xi(^) corresponds to the experimental macroscopic hnear response while Xaí^ij^^j^s) is the expe- 
rimental macroscopic nonhnear response. 

It is shown in ref. [2J], that for a field E{t) = Ecos{ujt) one gets : 



^{t) 
eo 



.(1) 



£:|Xi|cos(a;í-ái) + 3/4£;^ X3 cos{ujt ^ 6¡' > ) + 1 / AE"^ Xs cos{3Lüt - 6^^'') + 



(1)^ 



/3) 



?(3)^ 



(30) 



We now must identify the resuh of our model with above relations giving the experimental macroscopic polarisation. 
We start from Egs. í^Dll^ and average over 61 which yields : 



, , ^3 Ncorr{lJ'molec)^ ((cosé'i)-*) 9A^^, 

< MF >= — — :::i^ -E = 



ela\^Xi)\,_ 



a^{kBT)3 5(1-^2)2 keT 

We then average over the r's, as in Eq. [23l and then over Ncorr- With Eqs. Í2TII3Ü1 we obtain : 



(31) 



00 

/ 



[< MF"^ >] (1 - ¿2) . 

-^ ^ ^°" ^D^^^(wt)cos Swí + í-^^'^wr^-arctan^^wr) ^(hiT^dlnT 

4^1 + (3wt)2 L 



= (l/4)eo-B3 



(3) 



cos(3wí — ág ) 



(i/4)í; 



3e|(Axi)^ 
kRT 



Nc. 



\'H{ujTa)\ cos{3ujt + aTg{'H{ujTa))), 



(32) 



(3) . . ■ . (3) 

where the last equahty was obtained by replacing X3 by Bouchaud-Biroh's prediction Xs 



fcsT 



[Ncorr]av'H. {ujTa) , scc thc maiu articlc. Combining Egs. 1311(5^ one obtains 



11 



[NcorAav |H(wTa)|cos(3a;í + arg(H(wTa))) = 

- - /o\ r /n\ 

í?3 {lüt) cos Sojí + 5*3 {lüt) — arctan(3a;T) 



9[iVe. 



5(1 -¿2) 



^(In- 



v/1 + (3^t)2 



dln^ 



(33) 



.(3) 



.(3) 



As in the main article one defines the dimensionless nonlinear suscceptibility as X'í^^^ = x)im^BTl[eoar{¿^Xiy 
Writing X^l^ = [x'Í')n[ exp[-¿(5^^i] one obtains : 



and 



IX 



(3) I 

3,m I 



^(3) 
"3,m 



with S, 



(3) 

cos 



and S 



9[7V, 



corr\av 



5(1 -¿2) 



(Ai) \ , (o(3) V 



phase of X^^ = arctan 



(3) 



s 



(3) ' 
SIN 



(3) 

cos , 



s. 



(3) _ 



^(In' 



e(lni 



j-Q^ r j-oN 

I?3 (cjr) cos ^^3 (cjr) — arctan(3a;r) 



/Q\ r /-qn 

1^3 (c(jr) sin ^3 {lüt) — arctan(3cjr) 



\/l + (3c^t)2 



dln^ 



dhi'; 



(34) 



A similar calculation for Xg ^ yields : 



and 



IX 



(1) I _ 'j[-^*corrJaD 

5(1-^2) 



3,ml 



x(i) 

-'3,m 



('^COSJ + ('^S/Af j 



/5(1) ~ 

phase of X3 „j = arctan ^^J^ 



-"cos , 



with SfjQg = / ^(hiT) 



1)3 {ujt) cos ^3 (wt) — arctan(ajT) 



and S 



(1) 

S/JV 



e(hiT)- 



v/1 + (^t)2 

D3 ' {üjt) sin í'g ■^ {ojt) — arctan(a;T) 



v/1 + {^ry 



dlnT 



dln-; 



(35) 



Note that in the first equahty of Eq . [35l there is a factor 3 instead of 9 found in Eq. |34l This comes from Eg.lSOlwhere 
there is a factor 3/4 for the cubic term oscihating at w while it is only 1/4 for thc cubic term osciUating at 3lü. 



II. MORE ON THE FREQUENCY DEPENDENCE OF THE WEIGHT p = exp[ /e^//]. 

Ah above calculations have been made as if the hfetime t^x of the considered Asymmetric Double WeUs is infinite. 
As a supercooled hquid is ergodic above Tg, the heterogeneity of the dynamics imphes that t^x must be finite. This 
comes from the fact that a region of space relaxing faster than the average must become a region relaxing slower 
than the average, to restore ergodicity. We shaU assume, for simphcity, that any ADW is reshuffled with the same 
characteristic time t^x, whatever the value of t it had just before. 

After reshuShng, the glassy correlations are different from those estabhshed before. Thus, if one performs an 
average over time longer than t^x, a given molecule is no longer correlated to any other molecule. This is why, one 
expects any molecule to become effectively independent of aU other molecules in the hmit of large times t ^ Tex- 
Therefore one expects, at large times, the measured nonhnear dimensionless susceptibihties X3 to be dominated by 



the corresponding susceptibihties X. 



(fe) 

S.trivial 



of independent molecules undergoing Brownian rotational motion. Note 
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that Xg ¿¿^¿^; has been fuUy calculated in Ref. [2^ . 

Very few things are quantitatively estabhshed concerning the reshuffling phenomenon ; even the value of Q = Tex/Ta 
remains a subject of discussions 2Q\. Therefore a detailed description of its impact on nonhnear susceptibihties is not 
available at present. This is why, we heuristicahy add the nonhnear susceptibhties Xg ^ given by our ADWP model 

(muhiphed by the weight p) to Xg ^^^^^g^i (multiphed by the complementary weight [l—p)). As thc hmit of large times 
í ^ Tex corresponds to low frequencies fTex < 1, we physicaUy expect that p vanishes in this hmit. The simplest way 
to express this idea quantitatively is to state that the weight (1 — p) of the trivial response is given by the probabihty of 
a reshuffling event happening during one E oscihation period of ^tt/w. It is reasonnable to assume that the probabihty 
of the reshufñing events are given by a Poissonian distribution {1/Tex) exp {—t/Tex), and therefore : 

1 - P = / exp {-t/Tex) — which yields p = exp {-fex/f) where fex = I/tsx- (36) 

This is the weighting function that has been used in Figs. 3-4 of the main article. Of course it plays a role only for 
the range / < /q as one has p{f > /„) ~ 1 since fex/fa < 1- 

III. AVERAGING OVER 5. 

For simphcity we have presented in the main article the resuhs of our ADWP model obtained for a single value 
of the dimensionless asymmetry d. One can generahse the results by averaging over S, at the cost of additional 

(k) 

parameters. To investigate this question, we computed the values of X^ /^ for 100 values of S hnearly distributed in 

the [0;0.99] intervaL We then averaged the complex values of A¿ ^ by a weight w^S). For simphcity we have used 
either a flat distribution which is non zero only between Smin > and Smax < 1 ; or a "gaussian" distribution where 
w{S) = C X exp [—{S — (5i)^/(2 x (¿2)^)] . Here C is the proper normahsation constant taking into account that S is 
defined only on the [0; 1] intervaL Note that Si is close to, but not exactly equal to, the average of S ; and similarly S2 
is not exactly its standard deviation due to the fact that S is restricted to the [0; 1] intervaL 
Two interesting features are worth noting in this averaging procedure over S : 

(k) 

• First, the values of X^ ^ plotted in Figs. 3-4 of the main article can be recovered with distributed values of S. For 

example, Fig. S[Ü]below shows the values of Si and of ¿2 that have to be chosen to recover the values of X^ ^ plotted 
in Figs. 3-4, by using a gaussian distribution. One sees in Fig. S[S1 that the value S = 0.60, chosen in the main article 
to fit the experiments without averaging over S, corresponds to the hmiting case of a gaussian distribution with a very 
smah standard deviation. Beyond Si = 0.60, one cannot recover the curves for Ag j^ given in the Figs. 3-4 of the main 
article. 

• Second, the shape chosen for w{S) can strongly change the resulting X^ ^ values. To investigate this point, 
we have fixed the two first moments of (5, and chosen accordingly the parameters Si,S2,Smin and Smax- It is found 
that Xg j^ can be strongly different for a gaussian weight and for a flat weight distributions. This clearly shows the 
strong importance of the interference effects, evoked in the main article, between the nonhnear susceptibihties of 
the dynamical heterogeneities corresponding to different values of r. These interference effects are strong enough to 
yield, e.g., a change in the log-log slope of X¡¡ mif I fa > 1) as weh as a change in the values of A3 m{fa) by a factor 
significantly different from 1 (i.e., larger than 2, or smaher than 1/2). We emphasize that the changes of X^ ,„(/) ^^^ 
in most cases different from those observed on X^ mif)- This is the reason why it is not surprising that fitting the 

measured values of Ag and of A3 requires different values of [Ncorr]av, as in the main article. Indeed, it is very hkely 
that the extreme simphcity of our model cannot fuhy capture these comphcated interference effects. However, relaxing 
only this constraint that the values of [Ncorr]av should be the same when fitting X^ and when fitting X3 , we have 
shown, in the main article, that our ADWP model is able to reproduce the sahent features of the nonhnear experiments 
on glyceroL This is why we think that this model is reaüy relevant for showing what are the new informations about 
the glass transition that can be drawn from nonhnear experiments in supercooled hquids. 



13 



0,5 




0,0 0,1 



Fig. S 5: Values of ái, 5^ yielding, with a gaussian weight distribution, the same values of Xg ^ as those obtained in the main 
article with a single value of 5. The Une is a guide to the eyes. 



